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We set up a method for a recursive calculation of the effective potential which is applied to a cubic 
potential with imaginary coupling. The result is resummed using variational perturbation theory 
(VPT), yielding an exponentially fast convergence. 



o 

o 

(N 
(N 



> 

o 

o. 
o 

:^ 

Oh: 

rt ■ 



% 



PACS numbers: 02.30.Mv, 03.65.- 



I. INTRODUCTION 



Perturbation theory is the most commonly used 
technique for an approximate description of non-exactly 
solvable systems. However, most perturbation series 
are divergent and yield acceptable results only after 
resummation. In recent years, based on a variational 
approach due to Feynman and Kleinert [J, a systematic 
and uniformly convergent variational perturbation theory 
(VPT) has been developed H H S |B|. It permits the 
conversion of divergent weak-coupling into convergent 
strong-coupling expansions and has been applied suc- 
cessfully in quantum mechanics, quantum statistics, 
condensed matter physics, and the theory of critical 
phenomena. 

The convergence of VPT has been proved to be 
exponentially fast |3, Q, and this has been verified 
for the ground-state energy of different quantum me- 
chanical model systems. If the underlying potential is 
mirror-symmetric, one introduces a trial oscillator whose 
frequency Q is regarded as a variational parameter and 
whose influence is minimized according to the principle 
of minimal sensitivity p . In this way, the ground-state 
energy of the quartic anharmonic oscillator was analyzed 
up to very high orders in Refs. X-SJ- 

If the potential is not mirror-symmetric, the center 
of fluctuations no longer lies at the origin but at some 
nonzero place X. In VPT, this situation is accounted 
for by regarding the nonvanishing center of fluctuations 
X as a second variational parameter. An extreme 
example is a complete antisymmetric potential, such as 
V{x) = Ax^, which for real A does not correspond to 
a stable system. Interestingly, if the parameter A is 
chosen to be imaginary, so that there does not exist a 
classical system at all, the quantum-mechanical system 
turns out to be well-defined, and the spectrum of the 



Hamilton operator 



H 



2dx^ 



IX 



(1) 



is real and positive. This remarkable property of 
the non-Hermitian Hamilton operator, found in Refs. 
HIHIlllIHIl, can be attributed to the fact that it 
possesses a different symmetry: it is invariant under the 
combined application of the parity and the time-reversal 
operation. 

In this paper, we apply VPT to the T'T-symmetric 
Hamilton operator (Q. In a first naive approach, we 
ignore the necessary shift X of the center of fluctuations 
and resum the weak-coupling series of the ground-state 
energy for the anharmonic oscillator 



V{x) = —x^ + igx^ 



(2) 
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in the strong-coupling limit. In this limit, the potential 
reduces to the purely cubic potential of Q. It turns 
out that the VPT results approach the corresponding 
numerical value for the ground-state energy of (^ with 
increasing order, but the rate of convergence is not satis- 
factory. Afterwards, we allow for a nonvanishing center 
of fluctuations X by using the effective potential, whose 
calculation is accomplished by an efficient recursion 
scheme. This refined approach improves the convergence 
of the results drastically. 

In Section ^ we derive the weak-coupling series 
for the ground-state energy of by evaluating con- 
nected vacuum diagrams. In Section IIIII we show how 
this perturbation series can be obtained more efficie ntly 
by means of the Bender- Wu recursion method [lj|. 
In Section IIVI we resum the weak-coupling series for 
the ground-state energy of ((2Jl by applying VPT and 
examine the resulting convergence. In Section we 
determine the effective potential with the background 
method [l5|, llq from one-particle irreducible vacuum di- 
agrams. In Section IVll we set up new recursion relations 
for a more efficient calculation of the effective potential. 
In Section rVIII we finally treat the resulting expansion 
with VPT and examine the improved convergence. 



II. PERTURBATION THEORY 

The perturbation series for the ground-state energy of 
the anharmonic oscillator ^ can be calculated from con- 
nected vacuum diagrams. Up to the fourth order in the 
coupling constant g, the ground-state energy is given by 
the Feynman diagrams 
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with the propagator 



Gu{n,Tj) = — e 



-w Ti-r,- 



and the vertices 



--Jo '^- 



Evaluating the Feynman diagrams lj2l leads to the fol- 
lowing analytical expression for the ground-state energy: 



E = huj 



1 Ug^n 465/^2 



8w5 



32^10 



+ 0{g') 



(6) 



Since evaluating Feynman diagrams of higher orders is 
cumbersome, only low perturbation orders are feasible 
by this procedure. If we want to study higher orders, we 
better use the Bender- Wu recursion relations il^ ■ 



III. BENDER- WU RECURSION RELATIONS 

The Schrodinger eigenvalue equation for the anhar- 
monic oscillator (0), 

- Y^'"{^) + ( Ya;2 + tgxA i;{x) = Ei^ix) , (7) 



k 


1 


2 


3 


4 


5 


(■k 





11 

8 





465 
32 





k 


6 


7 


8 


9 


10 


(■k 


39709 
128 





19250805 
2048 





2944491879 
8192 



TABLE I: Expansion coefficients for the ground-state energy 
of the anharmonic oscillator (|5J up to the 10th order. 



The (pk (i) are expanded in powers of the rescaled coor- 
dinate x: 



fc+2 
m— 1 

(4) For the ground-state energy we make the ansatz 



(10) 



E^nw{- + Y,9''^k\ 



(11) 



(5) Inserting JHJ - (|ll|l into ((JJ), we obtain to first order 



„(i) 



.(1) 



0, c 
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-, ei-0. (12) 



For k > 2, we find the following recursion relation for the 
expansion coefficients in (I10|) : 
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(m + 2)(m+l) (fe) 



2to 



-m+2 



(13) 



- fc— 1 m+l 
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Jk) 



with Cm = for m > k + 2. The expansion coefficients 
of the ground-state energy follow from 

k~l 

(14) 



(fc) 



lE^'^ 



(0 Jk-l) 
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TablelHshows the coefficients e^ up to the 10th order. We 
observe that all odd orders vanish for symmetry reasons 
and that the first two even orders agree, indeed, with the 
earlier result ((HJ. 



is solved as follows: We write the wave function ipix) as 



'^i^) = [:^) <^^p 



'Y + '^(^) 



(8) 



with the abbreviation x = x^fujlh^ and expand the ex- 
ponent in powers of the dimensionless coupling constant 
g = g^/hJuJ^ by using 



^g'^Mx) 



(9) 



fe=i 



IV. RESUMMATION OF GROUND-STATE 
ENERGY 

In this section, we consider the strong-coupling limit 
of the perturbation series Hll() . Rescaling the coordinate 
according to a; ^ xg~^'^, the Schrodinger equation (jjj) 
becomes 

-Y& V'(x) + l^lg-'/'^'^' + ^x''^ ^ix) = (15) 



Expanding the wave function and the energy in powers 
of the couphng constant yields 



V'(x) = Mi) + r^'^'Mi) + 9~''"'Mi) + 



(16) 



and 



E 



huj g2/5 [bo + r^/'&i + .r^/'&2 + . . .) . (17) 



By considering (|15|l in the hmit g — > oo, we find that the 
leading strong-coupling coefficient 60 equals the ground- 
state energy associated with the Hamilton operator ^. 
A precise numerical value for this ground-state energy 
was given by CM. Bender |9lllq: 

&o = 0.762851773... . (18) 

The weak-coupling series (|11|) is of the form 



E'^^\a,uj) =huj 



N 

E 
fc=i 



ha\ 



^2k 



(19) 



with the abbreviation a 



'. Table U suggests that (O 
represents a divergent Borel series which is resummable 
by applying VPT H S S 0. To this end, an artificial 
parameter is introduced in the perturbation series, which 
is most easily obtained by Kleinert's square-root trick 



n^/T 



with 



-r!2 
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(21) 



Thus, one replaces the frequency m in the weak-coupling 
series ()19|l according to (|20|) and re-expands the result- 
ing expression in powers of a up to the order a^ . Af- 
terwards, the auxiliary parameter r is replaced accord- 
ing to (|21|l . The ground-state energy thus becomes de- 
pendent on the variational parameter fi: E^^'{a, uj) — )■ 
E'-^\a, uj, f2). The influence of fi is then optimized ac- 
cording to the principle of nrinimal sensitivity 6], i.e. one 
approximates the ground-state energy to iVth order by 



£;W = £;W(a, Lo, f]W) 



(22) 



where J7'^^ denotes that value of the variational pa- 
rameter for which E^^'{a, uj, fi) has an extremum or a 
turning point. 

Consider, as an example, the weak-coupling series 
((T^ to first order: 



^ ' ' 2 8^4 



(23) 



Inserting H20|l . re-expanding in a to first order, and tak- 
ing into account H21|l. we obtain 



Extremizing this and going to large coupling constants, 
we obtain the strong-coupling behavior of the variational 
parameter: 

rid) = ^^1/5 (^^(1) + ^7(1)^-2/5 ^ f^(l)^-4/5 _^)j ^(25) 

with the abbreviation a = g^ and the coefficients 
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Inserting the result H25|l. (|26|l into H24|l . we obtain the 
strong-coupling series (|17(l with the first-order coefficients 



(1) _ 5^22 
^° ~ 16 ' 
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The numerical value of the leading strong-coupling coef- 



ficient is b. 



(1) 



0.5799. Thus, to first order, the relative 



deviation of the result from the precise value lfTS|l is 
\b^i' - b, 



bo 



24% 



(28) 



Despite this relatively poor agreement, it turns out that 
the VPT results for 6q ' in higher orders converge to- 
wards the exact value (TSJ. In Refs. 01 3 it is proved that 
VPT in general yields approximations whose relative de- 
viation from the exact value vanishes exponentially. In 
our case we have 



h(^) 



bo 



(X exp (-CiV^/^) , 



(29) 



where the exponent 3/5 is determined by the structure 
of the strong-coupling series ((T7|l . 

In Fig. n tliG exponential convergence of our varia- 
tional results is shown up to the 20th order. Fitting 
the logarithm of the relative deviation to a straight line 
yields 



In 



60 



-0.96(ll)iV3/5 _ ^^83(44) (-3Q) 



In the following, we show how this exponential conver- 
gence is improved drastically by allowing for a shift of 
the center of fluctuations. 



V. DIAGRAMMATIC APPROACH TO 
EFFECTIVE POTENTIAL 

In the presence of a constant external current j, the 
quantum statistical partition function reads 



Z[j):^ jvx exp|-i^[a;](j)| , 



(31) 



where A[x]{j) is the Euchdean action: 



A[x]{j) = 



hl3 



dr 



-x\t)+V{x{t))-jx{t) 



(32) 



The free energy thus becomes a function of the external 
current: 



P 



(33) 



The path average, 

1 £ ^ \ r" d^ ^ ^ 

X = —— iVx I — - x(t) 
Z{j) J Jo h(3 ^ ' 



rhfS 



exp 



-j:AMj) 



(34) 



then follows from the first derivative of the free energy 
with respect to the external current: 



X 



dFjj) 
dj 



(35) 



Assuming that the last identity can be inverted to yield 
the current j as a function of the average X, one defines 
the effective potential T4ff (X) as the Legendre transform 
of the free energy with respect to the external current: 



V,siX)=F{j{X))+j{X)X. 



(36) 



Furthermore, the first derivative of the effective potential 
gives back the external current j: 



dX 



Kx). 



(37) 



Thus, the free energy F = F{j = 0) can be obtained by 
extremizing the effective potential. 



F = Kff (Xc) , 



with 



dV,s{X) 



dX 



= 0. 



(38) 



(39) 



X=Xe 



In the zero-temperature limit, extremizing the effective 
potential then yields the ground-state energy. 

The effective potential is usually not calculated by 
performing explicitly the Legendre transformation H3t)|) 
but by a diagrammatic technique derived via the so- 
called background method [I^ [l^ . There, the effective 
potential is expanded in powers of the Planck constant 
ft, and the expansion terms are one-particle irreducible 
vacuum diagrams. The result is 



V^s{X) = V{X) + -TrlnG-i + V^"''\X) , (40) 



where the trace-log term is given by the ground-state 
energy of a harmonic oscillator of X-dependent frequency 

Co=^W{X)- 



h „ T huj 



(41) 



The interaction term y(™*^(X) contains the sum of all 
one-particle irreducible vacuum diagrams. For the anhar- 
monic oscillator ^ , the relevant subset of the diagrams 
in © is 

y(int)/^^ = - lim kgri — fX + — 

+ 16 O } + ^^^'^ ■ (^^^ 

These one-particle irreducible vacuum diagrams are de- 
rived most easily by an efficient graphical recursion 
method [l3|- The frequency of the propagators is now 
given by: 



Cj = v/^2 4_ Q^gx . (43) 

By evaluating the diagrams (I42II we obtain 



V,s{X) = —X^ + igX^ + -^uj^ + 6igX (44) 
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The ground-state energy of the anharmonic oscillator Q 
is found by extremizing the effective potential 1)44(1 . To 
this end, we expand the extremal background according 
to 

X, = t{Xo + hXi + tex2 + h^Xs) + Oih*) . (45) 

Inserting (|45() into the vanishing first derivative of H44() 
and re-expanding in h, we obtain a system of equations 
which are solved by 



Xn = 



Xi = - 



_3g_ 

2tj3 ' 






Inserting H45|l . H46|l into (|44|l and re-expanding in h 
yields again the ground-state energy (|S|). 

In order to go to higher orders, we shall now de- 
velop a recursion relation for the effective potential. 



VI. RECURSIVE APPROACH TO EFFECTIVE 
POTENTIAL 

In the presence of a constant external current j, the 
Schrodinger eigenvalue equation for the anharmonic os- 
cillator |(2Jl reads 

fi,2 /(jj2 \ 

- y V'"(a;) + ( Y^^ + tgx^ ~ jx\ ij{x) = Ei;{x) . (47) 



Taking into account the Legendre identities (|36|) . (|37|l . 
Eq. ||T7|) becomes 

= V,s{X)-V:,s{X)x\4,{x). (48) 
If the coupling constant g vanishes, Eq. (|48|l is solved by 



V'(x) = AA exp ( Xi - y ) , (49) 



yeff(X) = nc.(i + | 



(50) 



where the path average has been rescaled by the oscilla- 
tor length: X ~ X^JuJJh. For a nonvanishing coupling 
constant g, we solve the differential equation (|48|l by the 
expansions 



ip{x) = TV cxp 
V,s{X) = hw 



Xx- — + <i){x) 



J + ^ + E^^-^^w 



/£=1 



(51) 
(52) 



For the correction to the wave function, 0(i), we make 
again the ansatz (jHl, H10|l . Thus, we obtain from H48|l for 
fc = 1: 



c^l^='- + 2^X^ 



.(1) _ 



Co = 






Ji) _ 



^iW = ^ + ^^^ 



-3 , (53) 
(54) 



For k > 2 one finds for m > 2 the following recursion re- 
lation for the expansion coefficients of the wave function 



Jk) _ im + 2){m + l) (k) X{m + 1) ,k) 



m 



k-l m+1 



-. k — i Tn-t-i 

— ^5:n(m + 2-n)cWc(t;L., (55) 



1 = 1 n=l 



(k) 

with Cm = for m > k + 2. For to = 1, we have 



(*-■) Q (fe) 



2Xc 



,(fe) 



^fe(^) 



+ E(4'"^4ncf-')c«) . (56) 

1=1 

The expansion coefficients of the effective potential follow 
from 



fe-i 



(57) 
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TABLE 11: Expansion coefficients for the effective potential 
of ij^J up to five loops. 



Using these results, the effective potential can be deter- 
mined recursively, yielding an expansion in the coupling 
constant g: 

VMX) = ^ + ^X^ + ^9('-^+X^) (58) 
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This result is in agreement with the expansion of l|44f) in 
powers of g and can be carried to higher orders without 
effort. The expansion coefficients for the fi-expansion 



N 



i4ff(x) = yn'y«(x) + o(n^+i) (59) 



;=o 



can then be obtained easily |1^] . Iterating the recursion 
relations H55|l - (|57|l up to the order g^, we obtain the 
effective potential up to five loops as shown in Tab. ^1 



VII. RESUMMATION OF EFFECTIVE 
POTENTIAL 

We now apply VPT to the loop expansion of the ef- 
fective potential (|59l) . Since the Planck constant h is 
now the expansion parameter rather than the coupling 
constant 5, Kleinert's square-root trick will be modified 
accordingly: 



fiVl + hr 



with 



LJ 



As an example, we consider again the first order: 
^X-^^gX-^- 



(60) 



(61) 



Kff^(^) = ^^' + ^9X^ + !^V^2 + 6,gX . (62) 



After substituting lo according to (|6()|l . re-expanding in 
?i, and taking into account (|61|l . we obtain 

v[§ {X, r!) = ^ a2 + igx^ + J yWTewx . (63) 



1-loop VPT 


0.742751023 


2-loop VPT 


0.764570478 


3-loop VPT 


0.758783545 


4-loop VPT 


0.762843684 


5-loop VPT 


0.762849959 


Numerical 


0.762851773 



TABLE III: Variational results for the ground-state e nerg y of 
Q compared to the numerical result 1181 of Refs. p.[lq|. 



In order to calculate an approximation for the ground- 
state energy, we now optimize in J7 and extreniize in X, 
yielding 



^v;S>,Y,o, 


x=x(i), n=n(i) 


= 0, 


(64) 


^KU'(.Y,0, 


x=x(i), n=o(i) 


= 0. 


(65) 



Equation H64|) is solved by 
Afterwards, we obtain from H65|l : 



(66) 



■iig 2y/6i^{XW)^/^ 



This equation allows us to determine the strong-coupling 
behavior of X : 



(68) 



V uj 
where the coefficients read 
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(69) 



Inserting the results (jHH), (EHl), EH into (E3) yields the 
strong-coupling behavior of the ground-state energy p7(l . 
with the new coefficients: 
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(70) 



The new numerical value of the leading strong-coupling 
coefficient is &q w 0.7428, which is in much better agree- 
ment with H18|) than the previous value of (|27|l : 

\^ll-M « 3% . (71) 
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FIG. 1: Convergence of the results for the strong-coupling co- 
efficient &o obtained from resummation of the weak-coupling 
series of the ground-state energy (circles). The resummation 
involving a variational path average X converges much faster 
(triangles) . The lines represent fits of the respective data to 
straight lines. 



Thus, the variational path average has led to a significant 
improvement of the first-order result. Table UTTI summa- 
rizes our results for bg ' up to the fifth order. Figure Q] 
shows the much faster exponential convergence up to the 
fifth order. A best fit of the data yields 



J.W 



In 



\h^ ' — h I 

^^^- ^ = -5.8(1.6)iV3/5 + 3.0(3.0) (72) 



6o 



(67) which is to be compared with H30() . 



VIII. CONCLUSION AND OUTLOOK 

We have developed a recursive technique to determine 
the effective potential, which is far more efficient than 
diagrammatic methods. In combination with VPT, this 
leads to a fast converging determination of the ground- 
state energy of quantum-mechanical systems with non- 
mirror symmetric potentials. It will be interesting to 
analyze in a similar way systems with a coordinate de- 
pendent mass term, where only a lowest order effective 
potential has been calculated so far |2(]j. Interesting fu- 
ture applications will address the effective potential of 0^ 
theories in 4 — e dimensions to obtain equations of state 
near to a critical point. A first attempt in this direction 
isRef. ISl. 
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